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The perturbative renormalization group for light-front QCD Hamiltonian produces 
a logarithmically rising interquark potential already in second order, when all gluons 
are neglected. There is a question if this approach produces also color van der Waals 
forces between heavy quarkonia and of what kind. This article shows that such forces 
do exist and estimates their strength, with the result that they are on the border of 
exclusion in naive approach, while more advanced calculation is possible in QCD. 
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In the constituent models used to describe hadrons, quarks can be assumed to 
interact via additive 2-body confining potential with the same color structure 
as if the interaction were mediated by one-gluon exchange, 



where F { ■ Fj denotes J2a=i FfF?, and Ff are the eight SU(3) color genera- 
tors for 2-th quark or antiquark, located at r-j. The potential models describe 
properties of single-hadron states quite well, but also predict power-law color 
van der Waals forces between colorless hadrons, which contradict experiments 
by being too strong. The color confining potentials also cause that the model 
Hamiltonians are unbounded from below for all states except color singlets, 
triplets and antitriplets |IJ. Potential models which try to avoid these draw- 
backs appear to be arbitrary and not systematically related to QCD. 

One of the attempts to derive a constituent picture of hadrons directly from 
QCD is by using the perturbative similarity renormalization group concept for 
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light-front Hamiltonian [§]. Perry [§], and Brisudova and Perry have shown 
that one can get confining potential already in second-order calculations based 
on this concept. In their calculation the model Hamiltonian is bounded from 
below for all color states of quark-antiquark pairs. However, their approach 
was not boost invariant and it was not clear how to describe two-meson states 
with clear control on their masses. Since then, the similarity renormalization 
group approach has advanced so that one can preserve boost invariance and 
required cluster properties by introducing effective particles [§, |J. 

We use similarity renormalization approach (see Ref. for description and 
references) to express the field-theoretical canonical Hamiltonian in terms of 
creation and annihilation operators for effective particles, which are unitarily 
equivalent to the canonical ones, viz. b\ = U\b CAn U x etc. The Hamiltonian H\ 
written using the new operators is band-diagonal, i.e. each term is multiplied 
by formfactor fx, which vanishes when energy changes by more than width 
A. We construct similarity transformation U x using second order perturbation 
theory. 

Following 0, we start from the canonical QCD Hamiltonian on the light front 
(LF), taking the advantage of the fact that in this formulation vacuum can 
be made trivial by means of a cutoff on longitudinal momenta. We regulate 
the Hamiltonian term by term in the expansion into products of creation 
and annihilation operators ||, introducing regulating factors r^si^,^ = 
exp k 12 /(iA 2 ) rs(x) for each operator in every term. Here x and fir 1 are 
the relative momenta of a particle with respect to other particles in a vertex 
(see ||). In gauge-field theory we need r$, the small-x regulator, and we choose 
rs(x) = Q(x — 5) as the simplest option. We use boost-invariant formfactor || 

h{Ml - Ml) = exp [-{Ml - Mlf/X 4 ] . (2) 

Here M 2 C = feiscreat Pj) 1S the free invariant mass of the particles created in 
the vertex, and Ml, similarly, for particles annihilated in vertex. 

In the evaluation of H\, we have to calculate the effective mass m x for quarks 
and antiquarks. Following [5|], we use the eigenvalue equation for a single quark 
state, derived in 2nd order for H\, at one arbitrary scale A = Ao to specify 
notation for the counterterm. 

™| 5 (P» = ^|«- ff ,™ ft ,| !(P , (3, 

Here M 2 is the free invariant mass of quark-gluon states. We express m\ Q , 
presumably infinite when the cutoff 5 is removed, in terms of the would-be 
perturbative eigenvalue m 2 . With hindsight, we set the eigenvalue mass rh to 
be equal to the constituent quark mass. Note that this condition is not meant 
to imply that there are free quarks, since it is imposed only in the formal 
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expansion in g, and ignores nonperturbative effects that prevent quarks from 
being free in this formulation. 



We use m\ as the starting point in evaluation of m\. Obtained quark mass 



m' 



+ 8m x diverges logarithmically when S goes to 0, 
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Here Cf is the SU(3) color factor for color triplet. We would obtain the same 
result for 8m\ in our approach if we used the coupling coherence condition 

10- 

Let us now consider the quark-antiquark interaction term in the renormal- 
ized Hamiltonian H\. For simplicity, we assume that all quarks have the same 
masses but different flavors (to avoid the insignificant discussion of antisym- 
metrization in the later treatment of the van der Waals forces). We will also 
use non relativistic (NR) approximation, which much simplifies further calcu- 
lations. This approximation is justified here when studying the bound states 
of heavy quarks because the formfactor f\ with A <C rh limits quark mo- 
menta to nonrelativistic values despite a large coupling constant f§. We use a 
three-vector qi2 defined through the requirement that 
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From this requirement one obtains x as function of qf 2 and q^ 2 , and one 
expands q\ 2 in the powers of |qi 2 |/m, obtaining in first order q\ 2 ~ (x — ~) 2m 
We should express the bare mass m in terms of m, but because the difference 
between m and m is of second order in g (see Eq. (4)), and the interaction 
term is already of second order, we can put m = m\ = rh here. 

When one applies the NR approximation to the effective potential f\V q q in 
Hamiltonian H\, one obtains the following result (cf. 0) 



A^NR(q) = -g 2 Fi2 ■ F 3 4exp 
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where q = qi 2 — q34. The NR potential is spin-diagonal. 
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, (6) 



The analysis of this interaction is made somewhat complicated by the external 
similarity factor fx, which depends not only on the momentum transfer q, but 
on |q_i2 1 and | Q34 1 separately. However, when A is much larger than momentum 
width of a wavefunction, fx can be replaced by 1 in the region that matters. 
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Following (|, |[], we neglect in our analysis of van der Waals forces all Fock 
sectors except QQ and QQQQ. This is justified for heavy quarks by the ex- 
pectation that gluons are lifted up in energy by certain gap condition, which 
is a non-perturbative effect. According to || fj, gluons should not appear 
explicitly in the resulting model. 

When gluons are neglected, there remains uncanceled (^-divergent, A-dependent 
term in effective quark (antiquark) mass term. This term cancels exactly 
in the matrix element |\I/2) against the divergence coming from QQ- 
interaction, when and |^ 2 ) are color singlets. For other color states, 
the Hamiltonian matrix elements are positive infinite, since the mass di- 
verges stronger than the potential of Eq. (6) 0. Therefore, the Hamiltonian 
is bounded from below for all states (as opposed to phenomenological additive 
color potential model where |\l/ 2 ) goes to — oo with growing distance be- 
tween quarks when both are in color-octet state, for all V{r) > oo [||]). 

This is the first important distinction between the renormalization group ap- 
proach to QCD and potential models, where the colored states would cause 
trouble. 

Note however that this cancellation does not occur outside the NR approxi- 
mation |J. One has to include transitions to qqg sector to obtain cancellation 
outside NR approximation. The second order perturbation calculations would 
then also lead to reduction of the confining part of effective potential, but non- 
perturbative effects could prevent complete cancellation |]X0|| , resulting only in 
weaker force. 

Color van der Waals forces between two mesons come from the mixing of 
octet-octet state (coupled to overall singlet) to singlet-singlet state, because 
the color dependent interaction polarizes mesons in color space as well as in 
position space. We choose states \a) = | I13I24) (both mesons in singlet state) 
and \(3) = | I14I23) , as the basis of the color subspace, where mesons are in 
color singlet state as a whole, where 1^- means that quark i and antiquark j 
are in color singlet. 

We will use static treatment to see if small-x divergences (which are constants) 
cancel in the eigenequation for QQQQ sector. For the meson-meson states, 



the interaction term Hj is a 2 x 2 matrix in color space |11 



Hi\a) = (§w„ - -up + ~u^\a) + (up - u q )\/3), (7a) 
Hi\P) = (u a - Ola) + (§ U/9 - \u a + |ii q )|/3), (7b) 



where 



U a = U 13 + U 24 , Uf3 = U l4 + U 2 3, U q = U l2 +U M , (8) 

and Uij = VijCij) = Vfa — rj) is the quark-(anti)quark potential. The mass 
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term is diagonal in |a), \/3). The divergences in the off diagonal (mixing) part 
cancel out because Vqq = V qq = —V q q. In the diagonal part, divergences in 
masses cancel with divergences in potential inside singlets (i.e. u a and Up, 
respectively) and divergences in the u q — Up and w q — u a cancel also. Thus, 
all small-x divergences cancel out and Eq. (7) has finite elements. Therefore, 
color van der Waals forces emerge in the approach presented here as they 
appear in color potential models. This is the second result of the effective 
particle approach in the approximation proposed in || [|. The third result is 
that we may then use the results on color van der Waals forces as known from 
potential models. 

To see what the van der Waals force looks like we have to look at the NR 
potential (6) behavior at large distances. We will use the spatial derivatives 
of the confining potential fl3j. The partial derivative with respect to r x for 
r z = fy = 0, r being the distance between Q and Q in quarkonium state, is 
given by the equation 

d a A 2 d 

T—V con{ (r x ) ~ 2V2tt— \og\r x /r \ + short range (9) 

or x 2ir m or x 

(the result for the r y direction is identical). Similar calculations of the shape of 
the potential in the r z direction give in the large distance limit the same func- 
tional form of the potential but twice weaker. Equation (9) contains unknown 
value of r , but from the fit to numerically obtained Fourier transformation 
one gets r of the order of m/A 2 . 

To simplify further calculations, we will use V con f averaged over angles (like 
in 0). Numerical calculations show that with growth of r the range of polar 
angle (9, where V C0Q f differs significantly from the value for O = ir/2, gets 
smaller and smaller. Therefore, for the large r we can use the analytical result 
for the i\l direction (9), getting the potential of the form 

u(r) = V Q log(r/r ), where V Q = 2a s \ 2 / [v^irm]. (10) 



For such logarithmic potential one obtains from the eigenvalues of 2 x 2 ma- 
trix in Eq. (7) and virial theorem QTJ] the following color van der Waals force 
between two heavy mesons: 

KMR) = h V ° {r R* \og(R/r 13 y (U) 

where R is the distance between the mesons and r 13 is the mean width of a 
meson. Note that the value of r does not enter the expression for the van der 
Waals force (11). 

The fit to the B meson spectra via heavy quark effective theory in || and the fit 
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to the cc and be spectra in |12| give a s X 2 /m around 1 GeV, i.e. Vq ~ 0.8 GeV. In 
color potential model (1) with V(r) = Volog(r/r ), one obtains from charmo- 
nium data V = 731 GeV [13|, [14). The question is then if such well established 
value of Vo produces acceptable strength of color van der Waals potential (11). 

We have done our calculations of the effective confining potential in the NR 
limit, which strictly speaking can be valid only for heavy quarkonia. To the 
author's best knowledge there currently exist no data on forces between such 
mesons. Therefore, we will use existing experimental limits on the color van der 
Waals forces between light hadrons, to make only order of magnitude guesses 
about the allowed strength of V . 

Between two (heavy) baryons one obtains the same type of force as between 
two (heavy) mesons but with different coefficient (e.g. Gavela et al. |15| gives 
1/2 for baryons instead of 3/54 for mesons, but they use a different method 
to calculate van der Waals force coming from potential (10)). Our NR poten- 
tial (6) does not apply to light quarks, but because van der Waals force is 
a low-energy interaction acting between slow-moving particles, and in con- 
stituent models NR approximation is commonly used, we will treat the NR 
equation (11) for the van der Waals potential as providing approximately 
right order of magnitude, disregarding any doubts that such procedure may 
no longer be adequate in the case of effective QCD for light quarks, where 
precise comparison with data is attempted. 



Data from Cavendish-type experiments |T6| give no real constraints on the van 
der Waals force resulting from logarithmic inter-quark potential, as opposed 
to the one from linear potential [I]]. The data from hadronic atoms ||16|| implies 
(3/54) V < 2 MeV i.e. V < 40 MeV for r 13 = 1 fm. We see that V ~ 0.8 GeV 
appears to be very large. Sawada JlT], finds in the behavior of P-wave 
amplitude of n-n and S-wave p-p scattering at low energy some signs of 
long range super-strong force oc 1/R n with n about 6, but no trace of forces 
oc l/[i? 4 log(-R)] is found. While the strength of his force is too large to be 
explained using a s that appears in Eq. (10), and remains to be explained, it 
is interesting that van der Waals force oc 1/R 6 comes from Coulombic part of 
the inter-quark potential, which cannot be easily canceled. 



It is clear that the calculation reported here can be considered only an ini- 
tial step towards derivation of color van der Waals forces in QCD. The non- 
cancellation of small-x divergences without NR approximation, and the lack 
of rotational symmetry of V q q, suggest that one needs to consider Fock sectors 
with additional gluons. Most prominently however, one should expect correc- 
tions from creation of quark-antiquark pairs ||19|| . These elements can reduce 



the van der Waals forces we obtained here and recover agreement with exper- 
imental findings. The most attractive feature of presented approach is that, 
contrary to the potential models, the effective particle picture is clearly open 
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to further study in QCD, and can incorporate all these effects starting from 
first principles. 
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